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Abstract 

We use wave packet mode quantization to compute the creation of massless scalar 
quantum particles in a colliding plane wave spacetime. The background spacetime rep- 
resents the collision of two gravitational shock waves followed by trailing gravitational 
radiation which focus into a Killing- Cauchy horizon. The use of wave packet modes 
simplifies the problem of mode propagation through the different spacetime regions 
which was previously studied with the use of monocromatic modes. It is found that 
the number of particles created in a given wave packet mode has a thermal spectrum 
with a temperature which is inversely proportional to the focusing time of the plane 
waves and which depends on the mode trajectory. 







1 Introduction 



Exact solutions representing the head on collision of two gravitational plane waves are some 
of the simplest exact dynamical spacetimes. They provide clear examples of highly non 
linear behavior in general relativity: when two plane waves collide the focusing effects of 
each exact plane wave lead to mutual fucusing. This is revealed by the formation of either a 
spacetime singularity or of a non-singular Killing- Cauchy horizon at the focusing points of 
the two waves |2|, [| . They may also be useful to provide local models for processes that 
may be taking place in our universe as a result of gravitational waves produced in black hole 
collisions |4], H, the decay of cosmological inhomogeneous singularities |J, or by travelling 
waves in strongly gravitating cosmic strings f?|, ||. 

Quantum effects in such dynamical spacetimes must surely be important and one ex- 
pects particle production and vacuum polarization when a quantum field is coupled to such 
background. Yurtsever || was the first to study field quantization on a colliding wave back- 
ground, he considered the Kahn and Penrose [p]] solution which may be interpreted as the 
collision of two impulsive plane waves. The solution which has curvature singularities at 
the focusing points of the plane waves allows for the definition of two physically meanigful 
vacuum states: an "in" vacuum state associated to the flat space before the collision of the 
two plane waves and an "out" vacuum state related to the flat spacetime regions behind 
the shock fronts also before the collision. The Bogoliubov coefficients relating the "in" and 
"out" creation and annihilation operators could be found only approximately in the long- 
wavelenght limit. In this approximation the spectrum of created particles is consistent with 
a thermal distribution. 

In a recent paper |9j we have considered the quantization of a massless scalar field in the 
background of the collision of two plane waves which form a non singular Killing- Cauchy 
horizon. This solution describes the collision of two gravitational shock waves followed 
by trailing gravitational radiation flnj. The interaction region of the two waves is locally 
isometric to a region inside the event horizon of a Schwarzschild black hole with the Killing- 



Cauchy horizon corresponding to the event horizon [fj, [11], [I2j. Two unambiguous and 
physically meaningful quantum vacuum states may be defined: an "in" vacuum associated 
to the positive frequency mode solutions in the flat region before the collision of the waves and 
an "out" vacuum related to the positive frequency modes defined through the two null vector 
fields in the Killing-Cauchy horizon. Such state, which is invariant under the symmetries 
associated to the horizon corresponds to the unique preferred vacuum state defined by Kay 
and Wald [O] in spacetimes with bifurcate Killing horizons. It was found that the "in 



n 



vacuum contains a number of "out" particles which is proportional to the inverse of the 
frequency. In the long-wavelength limit the spectrum is consistent with a thermal spectrum 
at a temperature which is inversely proportional to the focusing time of the gravitational 
plane waves, in agreement with Yurtsever's result ||. 

Our result however is exact: although the Klein-Gordon equation in the interaction 
region cannot be solved exactly the "in" modes become blueshifted towards the trailing 
points of the waves and can be propagated through this region by using the geometrical 
optics approximation. This is somewhat similar to the situation in a Schwarzschild black 
hole 0. 



In this paper we want to reconsider this problem in the light of wave packet mode quanti- 
zation instead of the monocromatic modes we have used in ref. PI. There are several reasons 
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that, we believe, justify this. One of the reasons is well known: the use of monocromatic 
plane wave modes leads to infinite expressions for the total number of particles created of 
a given frequency, whereas the number of particles created in a given wave packet mode is 
finite [111]. The second, and more important for us here, is that wave packet modes prop- 



agate in a simple way in a spacetime, such as ours, formed by the matching of different 
regions. Since the singular free coordinates of one of the regions differ from the next region, 
the propagation of modes which are extended through all space finds serious difficulties. In 
ref this could be done because of the blueshift of the modes in certain regions and the use 
there of the geometrical optics approximation which ammounts to ray propagation. Since 
wave packets are localized in space and the packet maximum follows a well defined path 
in the spacetime, a natural approximation may be taken by which the propagation of the 
packets through the different regions becomes a simpler problem. 

Another reason is that the wave packet formalism localize the fenomena of particle cre- 
ation: the one particle states defined with the wave packet modes have two labels, one gives 
information on the "energy" of the particle and the other on its "trajectory" (all within a cer- 
tain range of values). Finally, we know that our colliding wave spacetime can be maximally 
extended through the Killing-Cauchy horizon with the extended Schwarzschild spacetime 
|T0| , |T3J. This is possible if one of the transversal coordinates of the plane waves is made 
cyclic. The resulting spacetime represents the collision of two plane waves propagating in 
a cylindrical universe and the creation of a black hole of mass proportional to the strength 
(or focusing time) of the plane waves. In a forthcoming paper we want to consider particle 
creation in the extended spacetime. Using the results of the present paper the calculation 
will become somewhat similar to that of stimulated emission by black holes 16|| ; the use of 
wave packet modes has proved useful also in this case, see ref. jl7| . 

The plan of the paper is the following. In section 2 we introduce wave packet mode 
quantization. In section 3 we briefly review the geometrical properties of our colliding wave 
spacetime with special emphasis in the coordinates which are appropriate in the different 
regions. In section 4 we quantize a massless scalar field on the colliding wave background 
and propagate the "in" wave packet modes through the different regions. The advantage of 
wave packet modes is clearly seen in subsection 4.4 when the modes are propagated through 
the interaction region. The Bogoliubov coefficients relating packet creation and annihilation 
operators are derived in subsection 4.6, and the creation of particles is derived in subsection 
4.7 where we also compare our results with those of ref. |§. 



2 Wave packets 

Let us consider a complete orthonormal family, {/^(x)}, of complex solutions of the Klein- 
Gordon equation for a massless scalar field (f) (i.e. □</> = 0), which contain only positive 
frequencies with respect to a given timelike Killing vector d/dt, i.e. such that Cg/Q t fu(x) = 
—iujf^ix) with to > 0. Then it is possible to define a positive definite inner product between 
these solutions. The label u is continuous and stands for the energy (or frequency) and x 
stands for the spacetime coordinates. Since these solutions have a well defined value for the 
energy (or frequency) we can call them monocromatic modes. It is true that since the energy 
of the monocromatic modes is well defined their space localization is completely uncertain 
in concordance with the Heisenberg's uncertainty principle. 
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Now we want to work with a complete and orthonormal set of modes that are localized 
in space in some sense. To achieve this we can make an adequate superposition, within a 
small energy range, of continuous cu-labeled monocromatic modes in order to introduce an 
uncertainty in the energy and gain a certain information on spatial localization. We can do 



this as follows 14, 17, IS], define 



ye Jui 

where the new labels Co and n are restricted to verify that Co/e = j and ne/2-K = I are integers, 
e being a small and positive parameter. We will call this superposition a wave packet. Here n 
is a kind of Fourier label (the phase term exp(— inuS) has / periodes in the interval (Co, Co+e)), 
Co is the lower extremum of the integration interval and gives information on the energy of 
the wave packet. It can be shown easily that a set of discrete (Co, n)-labeled wave packets 
{fw,n{x)} is complete and orthonormal if the set of continuous ^-labeled monocromatic 
modes {fu(x)} is complete and orthonormal. 

It is worth noticing that a general set of wave packets, {fcj,n(x)}, given by ([l]), satisfy 
the following property, 

POO 

El/^)| 2 =/ doo\Ux)\\ (2) 
which follows from the equality, 

oo 

£ e ±in ^-^ = 65(00-00'). (3) 

l=— oo 

Note that the sums in (|[) and (H) are over the integer labels (j, I) which could also be used 
to label the modes, i.e. fj,i(x), instead of the labels (oo, n) which we will use throughout. 

We can now see in what sense the wave packets (|l|) are localized. Let us consider 
the generic wave packet (|l|) in terms of the modulus and phase of i.e. fu(x) = 

\Ux)\e ie ^\as: 

f , n (x) = ^ r +t dooe-™ +w ^ \f u (x)\. (4) 



If the interval of integration in ([]]) is small enough so that can be taken as approxi- 

matly constant (this will be true in general provided that has no singularities in the 

interval of integration), the 1/^(^)1 can be factorized out and we have an integral over the 
phase only. Now if the integrand's phase, i.e. Q u (x) = —noo + 9 UJ (x), oscillates rapidly over 
the range of integration (at least when noo is big enough, and this is true for I » 1) then 
the integral roughly vanishes except at the stationary phase points, that is when 



dGJx) 



doo 



0, (5) 



and we can use the stationary phase method to evaluate (01). Note that in the equation 



(j|) we have set oo = Co, after derivation, which is accurate provided the range of energy 
superposition is small enough. Equation (||) stablishes a relation between the labels (Co, n) 
of the wave packet and the spacetime coordinates x and so it represents a three dimensional 
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hypersurface. It determines the geometric locus of the spacetime points which give the main 
contribution to the integral (II]), this will be a spacetime region labeled by u and n and we 
will identify it as the wave packet trajectory, because outside this region the integral ([!]) 
roughly vanishes. It is worth noticing that if the monocromatic modes are labeled with 
some other continuous parameters, besides the energy, we can construct double or triple 
wave packets by superposition of monocromatic modes in a small range of these parameters. 
In this case we obtain surfaces or curves as the trajectories of the double or triple wave 
packet, respectively. 

Finally we can make the following remark about the uncertainty of the energy and 
position of a wave packet. When we construct a wave packet we take a superposition of 
Klein-Gordon solutions with a well defined energy u (the monocromatic modes) over a 
small interval (u, u> + e) , and so the uncertainty in the energy of the wave packet is Au ~ e. 
To know the uncertainty in the position of the wave packet let us assume the following 
simple form for the monocromatic modes fu(x) = e~ tuJt F w {x^) (where x\ j = 1,2,3 are 
space coordinates) which verify Cg/dt fu/{%) = —iojf^x). Then with the use of the stationary 
phase method we easily see that the wave packets (Q) are peaked arround values of time 
given by t = — 27r/e~ 1 , with width 2ne~ 1 (I is an integer) and so the time uncertainty is 
At ~ 2tt€~ 1 , which reflects Heisenberg's uncertainty principle Au; At ~ 2tt. 



3 Colliding plane waves geometry 

We will work in a spacetime that describes the head on collision of two linearly polarized 
gravitational plane waves propagating in the ^-direction. This spacetime has four regions 
(see Fig. 1): a flat region (or region IV) at the past, before the arrival of the waves, two 
plane wave regions (regions II and III) and an interaction region (region I) where the waves 
collide and interact nonlinearly. The geometry of these regions is given by the following four 
metrics, in coordinates which are adapted to the Killing vectors d x and d y of the spacetime 
by (see ref. H and references therein), 

t 2 a t t n • / \l2 7 7 1 — Sin (m + f ) , 2 

as-, = AL 1 L 2 1 + sin [u + v)\ dudv ; -ax — 

1 L 1 + sin (u + v) 



— [1 + sin (u + f)] 2 cos 2 (u — v)dy 2 , (6) 

(7) 



ds 2 n = ALiL 2 [1 + sin (u)] 2 dudv — y^--dx 

1 + sin (u) 



— [1 + sin (u)} cos 2 (u)dy 2 , 

ds 2 m = 4LiL 2 [1 + sin (v)] 2 dudv — 7— rcfe 2 

l + sm(u) ^ 

— [1 + sin (t>)] 2 cos 2 (v)dy 2 , 

ds 2 iy = 4LiL 2 dudv — dx 2 — dy 2 , (9) 

where u and v are two dimensionless null coordinates (v + u is a time coordinate and v — u a 
space coordinate) and L\, L 2 are two arbitrary positive lenght parameters, which represent 
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the inverse of the strength (focusing time) of the waves. The boundaries of these four regions 
are: {u = 0, v < 0} between regions IV and II, {v = 0, u < 0} between regions IV and III, 
{v = 0, < u < n/2} between regions II and I and {u = 0, < v < n/2} between regions 
III and I. 

At the surfaces u = n/2 and v = n/2 on regions II and III, respectively, the determinants 
of the respective metrics vanish, this marks the focusing points of the waves and a coordinate 
singularity. This singularity can be avoided with the use of appropriate coordinates (har- 
monic coordinates) in which the causal structure of these spacetime regions is well posed. 
In these coordinates the surfaces u = n/2 and v = n/2 become spacetime lines, see|| [UJ 
for details. 

Region I (in Fig. 1, this is the triangle bounded by the lines {v = 0, < u < vr/2}, 
{u = 0, < v < it/2} and u + t> = 7r/2)is locally isometric to a region of the interior of 
the Schwarzschild metric. This is easily seen with the coordinate transformation, 

t = x, r — M [1 + sin(w + v)] , 

if = 1 + y/M, 9 = 7r/2 - (u - v) 
where we have defined M = y/LxL^. The metric (§) becomes 

ds * = WL_{\ 1 dr i _ f _ dt 2 - r 2 [d6 2 + sin 2 ^^ 2 ) , 

which is the interior of the Schwarzshild metric. The surface u + v = tt/2 corresponds to 
the black hole event horizon. The boundary v = corresponds to r = M(l + cos#) and 
u = corresponds to r = M(l — cos^). These are the boundaries of the plane waves, 
these boundaries join at r = M (spacetime point of the collision) and also at the surface 
u + v = 7i/2ai6 = and 6 = it. This region of the Schwarzschild interior does not include 
the singularity r = and thus the interaction region has no curvature singularities. The 
above local isometry is not global however, the coordinates 9 and are cyclic in the black 
hole case but in the plane wave case, — oo < y < oo and — oo < v — u < oo. 

As in the the Schwarzschild case it is convenient to introduce a set of Kruskal-Szekeres like 
coordinates to describe the interaction region, because the (u, v, x, y) coordinates become 
singular at the horizon. These coordinates will play an important role in the quantization 
of the field. First we introduce dimensionless time and space coordinates (£, 77) 

£ = u + v, r) = v — u, (10) 

with the range < £ < n/2, —n/2 < i] < n/2 (we shall later see that in these coordinates 
the Klein-Gordon equation can be separated). Then we introduce a new time coordinate 
related to the dimensionless time coordinate £ by 

and a new set of null coordinates 

U = C-x, V = C + x. (12) 
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Note that the transversal coordinate x appears in the coordinate transformation because it 
behaves badly at the horizon. Finally, we define, 

U' = -2Mexp (- J^) <Q,V = -2Mexp < 0, (13) 

and the metric in the interaction region (|>D reads, 

= 2exp [(l- sinO/2] # , _ m2 )a 2 _ 

— (1 + sin£) 2 cos 2 ndy 2 , 

with 

U'V = 8M 2 T ^| 7 exp f fBLzl) , (14) 

(15) 

The curves £ = const, and x = const, are, respectively, hyperbolae and straight lines through 
the origin of coordinates (U' — V — 0), see Fig. 2. The Schwarzshild horizon (which is a 
Killing-Cauchy horizon for the spacetime) corresponds to the limit of the hyperbolae when 
£ — > 7r/2 i.e. the "roofs" V = or U' — 0. Notice that the problem with the transversal 
coordinate x at the horizon is that all the lines x = const, go through the origin of the 
(V, V) coordinates, so that all the range of x collapses into the point V = U' = 0, whereas 
the lines U' = and V = represent x = — oo and x = oo respectively. One should recall 
that we have not represented the coordinate x in our picture of the collision (Fig. 1) in 
which only the (u, v) coordinates are shown, x is a transversal coordinate perpendicular to 
the propagation and adapted to the Killing vector d x . 

To understand the global geometry of the spacetime one needs a tridimensional picture 
where the boundary surfaces between the different regions have to be written in terms of 
appropriate nonsingular coordinates adapted to each region. See [|, [HJ for details. 

It is worth noticing that most of the plane wave collisions produce true curvature singu- 
larities [Q] , but the spacetime described above is an example of a collision where the curvature 
singularity has been substituted by a Killing-Cauchy horizon (i.e. the surface U'V = in 
the Kruskal-Szekeres like coordinates). 



4 Wave packets in the collision of two gravitational 
plane waves 

4.1 Monocromatic mode quantization (summary) 

In ref. the quantization of a massless scalar field was considered in the spacetime 
described above representing the head on collision of two gravitational plane waves. Let us 
now briefly summarize this quantization sheme. We introduce a massless scalar field on 
the colliding spacetime background, which satisfies the Klein-Gordon equation, 

□0 = 0, (16) 
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where □</> = (—g) l ^ 2 {{—g) l ^ 2 g^ v 4 > .y) > an d fl 1 is the determinant of the metric. It has the 
generic plane wave solution 

4> = -=/(u, ^6^+^, (17) 

where the labels k x and fc^ are two separation constants which are interpreted physically as 
the momenta in the directions given by the Killing vectors d x and d y , which generate the 
plane symmetry of the whole spacetime. Introducing fllTD in ( f[6f) one obtains, 



{Vfg 



VFG 




f = 0, 



where the coefficients F, G and e , come from the adoption of a generic metric for all 
colliding plane wave spacetime, adapted to d x and d y , i.e. 

ds 2 = e~ Niu ' v) dudv - F 2 {u, v)dx 2 - G 2 {u, v)dy 2 . (19) 

One defines an initial vacuum state constructed with a complete orthonormal set of 



solutions of the Klein-Gordon equation (16) which are defined in the flat IV region, before 
the arrival of the plane waves, to be of positive frequency with respect to the timelike Killing 
vector d u+v ; these are the "in" monocromatic modes. Then one propagates these modes 
throughout the spacetime up to the horizon of region I (solving the appropriate boundary 
conditions imposed by the different classes of solutions for the Klein-Gordon equation, (|IBD , 
in the different regions of the spacetime). It is possible to define another "natural" vacuum 
state on the horizon of region I from a complete and orthonormal set of solutions of the 
Klein-Gordon equation (|16|), the "out" monocromatic modes, which are positive frequency 
solutions with respect to the vectors d/dU' and d/dV, which are two null Killing fields over 
the horizon. 

Comparing the propagated "in" monocromatic modes and the "out" monocromatic 
modes on the horizon via a Bogoliubov transformation, lead us to show that there is spon- 
taneous creation of particles in this spacetime with a spectrum of "out" particles given by 
the formula IP], 



•uuiin :..\ _ (2M) 3 <5 m> _„ f f f JU \C V 



(0, in|iV° ut |0, in) = ^ '- , m '~ m . / dk„ / dk (20) 

Here is the usual number operator for "out" particles, cj_, I' and m! are labels of the 

"out" modes (u;_ is the energy label), k and m are dimensionless labels of the "in" modes 
and (2M) 3 /7r|CV| 2 is a geometric factor whose coefficient, |CV|, depends on I', m', k_ and k. 
This spectrum is inversely proportional to the inverse of the energy of the "out" particles 
which are produced and it is consistent, in the long wavelenght limit, i.e. 8ttMuj^ « 1, 
with a thermal espectrum with temperature T = (87rM) _1 . The temperature is inversely 
proportional to the focusing time of the plane waves, given by the parameter M = \JL\Li-, 
and (4M) _1 is also the surface gravity of the horizon. Note that the spectrum contains 
a logarithmic divergence, this divergence appears because the spectrum which has units 
of (length) 3 , given by the factor (2M) 3 , describes the total number of particles created in 
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the whole spacetime volume. As it is well known [^3j this is characteristic of the use of 
continuous labeled modes and can be avoided using wave packets. 

The propagation of the "in" monocromatic modes throughout the spacetime is a difficult 
task. This is due, essentially, to the fact that the monocromatic modes are defined in the 
whole spacetime and there is not a single singularity free coordinate chart for all the four 
spacetime regions. Basically the problem is that the matching has to be done through all 
the spacetime points of the boundaries. Due to their spatial localization wave packet modes 
will be more easily matched. The matching will be done approximately on a single boundary 
point for each wave packet. 

In what follows we introduce the wave packet formalism in the colliding spacetime back- 
ground. As the first step we construct a complete set of "in" wave packets, from a superpo- 
sition of positive frequency monocromatic "in" modes, in the four spacetime regions. Then 
we will propagate these "in" wave packets throughout the spacetime up to the horizon, with 
appropriate matching conditions. These matching conditions will ensure that the trajectory 
of the "in" wave packet will be smoothly connected through the boundaries between the 
different regions. Next we will contruct a complete set of "out" wave packets on the surface 
of the horizon. Finally, we will relate the "in" and "out" wave packets via a Bogoliubov 
transformation and compute creation of particles in this formalism. 

4.2 Flat region (region IV) 

In the flat region the complete set of "in" -modes is, 



Uk x k v kAv>,v,x,y) = —j= e + v % (21) 

<2k^(2n) 3 



where the labels k x , k y and fc_ are independent separation constants for the Klein-Gordon 
equation fll6|) and where u', v' are two dimensional null coordinates related to the dimen- 
sionless null coordinates u, v by: v' = 2L2V, u' = 2L\U. The label k + is determined by the 
relation, 

Ak + k_ = k 2 x + k\. (22) 

It was shown in ref || that these modes are well normalized on the hypersurface {u = 0, v < 
0} U {v = 0, u < 0}. 

The labels k x and fc_ are continuous but k y is discrete if we take a cyclic spacetime in 
the ^/-direction (this is not necessary but it is convenient if we want to maximally extend the 
collidig wave spacetime later on). We identify k y with m/M where m is an integer. Our aim 
now is to discretize the continuous labels k x and k_ by constructing a double wave packet 
as, 

mn 1 rkx+t rk^+5 „-iqk--ink x 

ufl - k (u, v, x,y) = -= I dk x I dk. = e~ ik - v (23) 

y/eSJk. Jk. pkl{2nf 

Note that since we are only integrating over positive frequency monocromatic modes the 
vacuum associated with the wave packets is the same as that defined by the monocromatic 
modes (^TJ). 



8 



The integrand's phase is 



/£2 _|_ £2\ 

Q( 1v ^ = —nk x — qk_ — k_v' — ( ~~^r — ~ ) u' + k x x + k y y, 



(24) 



where the relation (|22|) has been used. The trajectory of the double wave packet, in the 
sense given by the stationary phase method, i.e. dQ^' /dk x = jdk- = 0, is given by, 



x^ = n + rJ-LiU, 
k^ 



2~U ' 



(25) 



(26) 



which define a null geodesic in the flat region. That is, this double wave packet moves on a 
null trajectory. 



4.3 Plane wave region (region II) 

The complete set of solutions of the Klein-Gordon equation (^) in region II is easily found 
(see ref. [[J) and is given by 

UkXk-( u , v , x iy) = I 1 3 — — ex P \ -hr (f( u ) k l + 9(u)kf) - ik-v' + ik x x + ik y y\ , 



^/2fc_(2vr) 



(27) 

where the two fuctions f(u) and g(u) are 

(l + sinw) 2 15 15 

j(u) = (9 — smu) H cosm u — 12, g(u) = tan it. (28) 

2 cos u 2 2 

The labels k x , k y and have the same meaning that in the flat region IV because this 
expression for the solutions of the Klein-Gordon equation in region II match smootly (i.e. 
in a continuous and diferentiable way) with the respective solutions ([H]) on the boundary 
between regions II and IV, i.e. {u = 0, v < 0}. 

In analogy with the previous case, we now construct a double wave packet 

(I « 1 rkx+e r-k-+S e -iqk--ink x 



k x ,n,k y ,k.,g - y/eSJk. J~k- ,/2A;_(27r) 3 

x — — exp f -i^r~ (f(u)k% + g(u)kT\ - ik-v' + ik x x + ik y y) . (29) 

pos ?/ \ ?.k: \ y > I 



1 I .L x 
exp -i— 

cosm y 2k 

The integrand's phase is given by. 



e(n) = _ nkx _ qk _ _ k _ v > _ (f(u)kl + g{u)k 2 y ) ^- + k x x + k y y, (30) 

and the trajectory of the double wave packet, i.e. dQ^/dk x = / dk_ = 0, is 

x^ = n + ^-L l f(u), (31) 
k_ 
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/(II) 



-q+(f{u)k 2 x + g(u)k 



2k 2 _ 



(32) 



which represents a null geodesic travelling throughout the plane wave region II. Note that 
this trajectory matches smoothly (i.e. in a continuous and diferentiable way) with the 
trajectory of the wave packet in region IV, i.e. (|25|), (|26|), on the boundary {u = 0, v < 0}. 
Since /(0) = g(0) = and df/du(0) = dg/du(0) = the following matching conditions are 
satisfied: 



x^ w \u = 


0) =x (n) (u = 


= 0); 




= 0) 




dx^ 


dxW 




dv'W 




dv'W 


du 


u =o du 


> 

u=0 


du 


u=0 


du 



u 



0), 



u=0 

This is not surprising because the monocromatic modes which we have used to define the 
double wave packets in the two regions, match smoothly on the boundary {u = 0, v < 0}. 

Before going further, we can extract some information on the meaning of the new discrete 
labels n and q. First, note that at the boundary between regions II and IV, {u = 0, v < 0}, 
we have from (g|), @ and (^TJ), (^2j), 



x =x (w) (Q) 



(ii) 



(0) 



x 



/(ii) 



(0) 



(0) = ~q, 



which means that a wave packet labeled by n and q has a peak on the null geodesic which 
crosses the flat region into the plane wave region at the spacetime coordinates Xq = n and 
v' = —q. Note that for wave packets in region III we can repeat the same discussion. 



4.4 Interaction region (region I) 

In this region mode propagation is more difficult; instead of relaying on the calculations 
given in ref. ||, we will start the discussion from the beginning. First, let us consider that 
the Klein-Gordon equation (|16D in this region can be separed by taking 

<f>(C,V,x,y) = e ik ^ +ik ^i; akx (C)^ aky (v), (33) 

and reduced to equations for ip a k x (0 anc ^ Vakyiv)- The coordinates £, r\ are related to the 
usual null coordinates u, v, by ([IT]), and the two new equations read, 



ttt - (tanO^ + ( a + f ^o^" ^ = °' (34) 



/ A; 2 l \ 

(p >m - (tmrf)(p jr ,+ [a- — =- \¥ = 0, (35) 

where a is a dimensionless separation constant and k x , k y are the same labels as in regions 
IV or II. We use the notation a = 2Ma, therefore k x , k y are dimensionless parameters. 

These differential equations have singular points at £ = tt/2, i.e. on the horizon of region 
I, and at r\ — ±7r/2 respectively. To avoid them we can perform the following change of 
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variables: <i£* = M(l + sin £) 2 cos x ^<i^, drj* = Mcos l rjdr]. With appropriate integration 
constants, £* reduces to (|TTD and 77* to: 



* = Mln 



1 + sin 77 

cos 77 



We can still introduce a new function 7^(0 instead of i/) a k x {£) by, 7 
and we arrive at the following equations for J a k x (£) an d <p a k y {v), 



fc 2 + 



cos 2 £ 



M 2 (l + sin0 4 



a; + 



2 COS 2 !] 

"^ + lvf 2 - a 



2sin£ 
1 + sin£ 



<^ = 0, 



7 



0. 



(36) 



(l/2)(l + sin0^ 



(37) 



(3? 



which have no singular points in the regions of interest. 

Recall that we are looking for a solution of the scalar field in the interaction region, 
restricted to satisfy certain boundary conditions imposed by the wave packets travelling 
from regions II (or III) into region I on the boundary {v = 0, < u < 7r/2} (or {u = 0, 
< v < 7r/2}) Thus the general solution for (ft in terms of the new functions j a k x (0-> ¥ak y (v) 
is given by, 

jk x x+ik y y \ f~i Tafc;c(£) 



<K£, v, x ,y) 



(39) 



1 + sin £ 

where the coefficients C a depend on a and the separation constants used to label the 
monocromatic modes in region II, i.e. k x , k y and /c_. Now we will try to obtain all possible 
information on the coefficients C a in ([39]). We know that expression (139|) gives the general 
solution (once k x and k y are fixed) for a massless scalar field in the interaction region, and 
the coefficients C a have to be such that the appropriate boundary conditions are satisfied. 
If we want to match the monocromatic modes defined in region I with those of region II (|27| ) 
(or III) we have to perform a sum over all possible values of a in (|39|) , this was one of the 
main difficulties in ref. 0. In the following discussion we will show that when we match 
the spatially localized wave packets travelling from region II (or III) into region I (^) with 
wave packets defined using the general solution (|3"9"D, there is only one coefficient C a which 
carries the main contribution. That is, we will find a single value of a for which the infinite 
linear combination (^) can be approximated by a single term. Furthermore, we will be able 
to find the value of the coefficient C a for that particular value of a. 

We will proceed as follows. First, we will find an approximate solution of equations (|37|) 
and (^) for 7^(0 and f a k y (v) respectively, near the boundary of region I with region II 
which will lead to an approximate solution of (|39D near this boundary. With this solution we 
will construct a wave packet in analogy with (|23| ) or (^9]) and then we will find the particular 
value of a and the phase of the coefficient C a which allows the matching of this wave packet 
with the wave packet (^) travelling from region II. Next, we will find an approximate 
expression for the general solution ( |3~9"D near the horizon of region I, i.e. the surface £ = tt/2, 
with the particular values of a and the phase of C a calculated before. Such solution becomes 
exact on the surface of the horizon and it will be used to construct a wave packet there. 

Let us start with the approximate solution of (j39|) near the boundary between regions II 
and I. First note that the solution <fi of the Klein-Gordon equation in the plane wave regions 
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(regions II or III) takes an exact WKB form, i.e. it can be written as <fr — CexpiS, where 
C and S are two real functions of the spacetime coordinates. This is directly related to 
the fact that the geometrical optics approximation is exact in the single plane wave regions, 
i.e. the rays of the Klein-Gordon solutions (the lines perpendicular to the constant phase 
surfaces) follow null geodesies. Thus we can expect that an approximate solution of the 
massless scalar field in the interaction region close to the boundaries with the plane wave 
regions can be obtained with the WKB method, this is physically related to the fact that 
near the boundaries the colliding plane waves superpose linearly. In fact from equations (]3Tf) 
and (|38|) we can see that they admit WKB solutions given by, see for example [p0~ |, 



7 wkb(D = -£= exp (±i f Q(C)di 

Jo e*) V 4* 



(40) 



¥wkb(v*) = i = ex P ( ±« / Q(V*) drf 
IQ{rf) \ J v* 



where Q(£*) and R(rj*) are, 



(41) 



Q 2 (D = k 2 x + 



cos 2 £ 



M 2 (l + sin£)' 



a + 



2sin£ N 
1 + sin£ 



t-,2/ *\ t 2 cos 2 ?? 



(42) 



(43) 



These WKB solutions can be used provided Q(£*) and R(rj*) do not vanish and they become 
accurate solutions when Q(£*) an d R{v*) change slowly, i.e. when \dQ(£*)/d£*\ « Q 2 {C) 
and \dR(rf) / dr)*\ « R 2 (r]*). Note that equation ( |3~7l) admits a WKB approximation of this 
kind throughout region I. Furthermore this approximation becomes asymptotically exact on 
the surface of the horizon of region I, i.e. £ = n/2, and it is a really good approximation 
near the boundaries between regions I and II (or I and III) provided the separation constant 
a is large enough: in fact, we will show that to ensure the correct matching between wave 
packets travelling from region II into I (or from region III into I) with a single term in the 
infinite sum fl39|) , it is necessary that the terms containing a in (pE2|) and (H), which also 
contain factors cos 2 £ or cos 2 ??, be dominant. Recall that when we aproach the horizon, i.e. 
the surface £ = n/2, on the boundary between regions I and II (or I and III), then both 
cos 2 £ and cos 2 r/ go rapidly to zero. This means that a must be big enough to compense the 
decreasing behavior of cos 2 £ and cos 2 ?]. Equation fl38| ) on the other hand admits a WKB 
approximation only near the boundary between regions I and II (or I and III) under the 
assumption of large a (this is because of the minus term in (f|3)). Fortunately, we know 
its exact solution when it is written in the form (|35|), this solution is given in terms of the 
associated Legendre polynomials and that will be the solution we will take in the regions 
where the WKB approximation is not valid. 

Let us now evaluate the form of 7^(0 and <fak y {v) near the boundary between regions 
I and II, i.e. {£ = -rj, < f < vr/2}. Although Q(«f) and R(r)*), as given by (g^) and 
(0), are the complete terms that appear in the WKB formulae, if we assume that a is large 
only the first two terms in powers of a -1 are relevant. With this expansion, performing the 
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integrals in (|40|) and ([y]) and choosing the minus sign in the exponent of ( f40"D and the plus 
sign in the exponent of fl4"I|), we have: 



_^ikxX-\-ikyy 



c <* i=k ~r exp ^ ~ ~ (^ - 6)] 

„ V cos 4 cos 77 L 



< (*S [/(C) - /(Co)] + fofo) - ^(770)]) + O (a- 3 / 2 ) |> , 



2v^ 



(44) 



where the functions / and g are given by (|28 



Next we construct a double wave packet as 



■.00 



,(£> »7, 2/) 



fcl+£ ^ x . [ k ~ +5 dk_e- iqk — ink ^ {l) 

Jk- 



k x ,n,k y ,k— ,q v 3 ' " ' " ' \J Jk 

and write the coefficient C a as and phase by 

C a = |C a | e l ® a . 

The integrand's phase in ([45"D is then, 

9® = — gA;_ — nA^ + fc^ + k y y + 9 a — 2(u — u )\^a 



n!k v kAt,V,x,y), (45) 



(46) 



M 2 



(kl [/(O - /(Co)] + A£ fofo) - g( V o)}) ^= + O ( 



-3/2 



(47) 



Now provides us with a natural way of matching the wave packets in region II and the 
wave packets in region I. Let us assume that the matching is good enough if the trajectories 
of the wave packets in regions II and I are joined in a continuous and diferentiable way on 
the point (u = uq, v = 0) of the boundary between these regions. This means that 9 (I) must 
satisfy the three constraints, 

(48) 
(49) 
(50) 



0(1) 


(u=uo, v=0) 


©(ID 


(« 


=uo,v- 


=0)' 






00(H) 






du 


(u=uo, v=0) 


du 




(u=UQ, 


v=0) 


39(1) 




00(H) 






dv 


(u=uq,v=Q) 


dv 




(u=UQ 


v=0) 



Using (|47| ) and fl30"D, these constraints give, respectively, 



9 a 



2^/a + (k 2 x f(u ) - k 2 y g(u ) 



Kf(uo) + k y g(u ) 
M 2 



2k. 



2y^ 



Kf( u o) + k y g(u ) 



2k. 



(k 2 J(u ) + k 2 y 
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M 2 
2JE 



2Lok„ 



(51) 

(52) 
(53) 



where f = df / du, g = dg / du are, using 



sinx^ 4 



f(x) = - =-^-, g(x) = — r -. (54) 



COS 2 X " COS 2 X 



Note that the constraint (|5T| ) gives the value of the phase 9 a for the coefficient C a in (p9[) , and 
the constraints (|52|) and (|53| ) define the value of a. These two constraints are compatible, 
this is not surprising because we have already noticed that the WKB approximation is the 
natural approximation near those boundaries. The value of a is, 



(klf(u )+k 2 y g(u ))-^-, (55) 



where we assume an expansion in even powers of cosmo and we recall that f(uo) and g(uo) 
are of the order (cosu )~ 2 . Then a goes like (cosu )~ 4 , this justifies our assumption that 
ctcos 2 £ » 1 or a cos 2 // » 1 because both terms go like (cosm )~ 2 in the matching point 
(« = Uq, v — 0). Note that (cosm )~ 2 increases rapidly when u differs from zero, but even 
when Uq = (|55"D with the use of (|22l) gives \fa ~ k + Li and the wave packet (]45f) still 
matches with (|29|). 

It is important to note that this matching fixes the phase of the wave packet but does 
not ensure its normalization. In fact, the normalization condition fixes the modulus of the 
coefficient C a , but it is convenient to postpone the calculation of this term until we have the 
form of the wave packet close to the horizon of region I, where the wave packet normalization 
condition is much simpler. 

Our next step is to identify the integrand's phase defining the wave packet (^) near the 
horizon £ = ir/2. To do this we can still use the WKB approximation fl40D of equation (|37|) 
near the horizon, this means that we can take asymptoticaly, 

7wKB (r) = ^e-^\ (56) 

\kx\ 

where we have choosen the minus sign in the phase in order to be consistent with (|44|). From 
(P^) the scalar field <f> is then given by 

( p-ilMCr-z;) ■ b > n 

0® &wy) = - 1 =T. c « Vc*, (v) >< ( 57 ) 

Vl^l a { e ~ llk ^ +x) ; k x < 0, 

where £* - x = U, C + x = V and U' = -2Mexp(-U/2M), V = -2M exp(-V /2M), 
following equations (0), Notice that when k x > 0, the scalar field reaches the "roof 

V = of the horizon (strictly speaking, the rays of 0, i.e. the lines normal to the constant 
phase surfaces of 0, reach the "roof V = 0), and when k x < 0, they reach the "roof 
U' = of the horizon. This asymptotic solution has been obtained because near to the 
horizon £ = 7r/2 the dominant term in fl4*2"|) is k x . However, it is not possible to obtain such 
an asymptotic solution for equation (0), because of the minus term in (|43|). Fortunately, 
as we have said, we can go back to the untransformed equation ( |35"|) which is the equation 
for the associated Legendre polynomials P a \~ (sin r\) . Note that in the cyclic case Mk y = m 
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is an integer and we can take a — 1(1 + 1) for I integer, then e lkyV P a k y (sinn) is the spherical 
harmonic Y™ (y/M, n/2 — n). The matching of the packets requires that a is given by (|5"5]). 
therefore we can take I = l a ~ y/a, where l a is the closest integer to yfa for large a. Now 
using the asymptotic form ([57]) of <fi near the horizon we can construct a double packet as 
in fl45]) . The integrand's phase defining the wave packet is given by, 

eW = - q k_ - nk x + k y y - k x (-x±C) + a , (58) 

where the upper sign in fl58[) stands for k x > and the lower sign stands for k x < 0, and 



9 a is the phase for the coefficient C a given in (|5"T|). As we have said the coefficients C a are 
restricted to satisfy a normalization condition, to make sure that the scalar field (f> is well 
normalized on the horizon. The appropriate inner product of any two solutions 0i, <p2 of the 
Klein-Gordon equation is given on the horizon by, see |J: 

r <-> 



J 2 ) = -i4M 



+ 



cos rj di] dy 

(01<9[/'</>2 



(7'=0 



(59) 



v=o 



Since, as we have just seen, 

C t (l,k a 



ky^ k — ) = oc e a 5 (/ / )) 



where l(k x ,k y ,k. 
(|59T) leads to, 



Z a and ^q, is given by (0), the normalization condition with respect to 

-6{l-l a ). 



40 o 



Ci (/, k x , ky, k_ 



(60) 



2M V 4tt£l 

Finally with the asymptotic value (|58| ) for the phase near the horizon (^ = 7r/2) we can 
use dQ^ /dk x = dQ^ /dk_ = to obtain the trajectory of the wave packet near the horizon: 



Q 



(k 2 x f(u ) +k 2 y g(uo) 



2k 



2 ■ 



X 



(I) 



A. i 



(61) 



(62) 



Notice that n + (k x /k^)Li f(uo) = x^^O) is the value of the coordinate a: when the wave 
packet trajectory coming from region II crosses the boundary v = into region I, as one can 
see from (|31|), then ( |62"D reads, 



x 



(i) 



x 



(I) 



(o)±r. 



(63) 



This is the equation for a null geodesic in a region close to the horizon as it can be seen by 
writing it in the form x 



(i) 



x 



(0) =p ln(cos£), where x^ '(0) is a constant for a fixed wave 
packet ((|Tl~D has been used). It is worth noticing that this is the equation for null geodesies 
close to the horizon as one can see in the Appendix of ||. 
Equation (^) can be written, using (0) and ([13]), as 

V = -2Me xm ^ 4M e 

cW(0)/4M 



U' Q ; k x > 0, 



V 



-2Me~ 



K ]k x <0. 



(64) 
(65) 



That is, the trajectories of the "in" wave packets, in the region near the horizon, are straight 
lines of U' = U' Q = const for k x > and straight lines of V = Vq = const for k x < 0. 



15 



4.5 "Out" wave packets 



Since at the horizon the fields d/dll' and d/dV become two null Killing vector fields we 
can define a new complete set of "out" modes on region I. We define solutions with positive 
frequency, uj + , with respect to the vector d/dU' on the V — "roof of the horizon and 
positive frequency, u;_, with respect to the vector d/dV on the U' = "roof of the horizon 
H [H| |21|. These modes, close to the horizon, are given by, 

uTvm>{U\V'^y) = e"** +u '~ l " ~ v ' , (66) 

2MJ(2tt)2cj_ 



where the labels m! and V are the three independents separation constants of the Klein- 
Gordon equation (|1(J) using the asymptotic metric on the horizon of region I, and the label 
u) + is given by, 

IQM 2 uj + uj_ =l'(l' + 1). (67) 

If we restrict, now, to the cyclic case the labels /', m' are both discrete but u_ is continuous. 
We can transform into a discrete label by constructing a single wave packet as follows, 

«SWM",1,»> = 4, f^.y'W'/ 2 -^^-^. (68) 

vV 2M v /(2vr)2^_ 

The integrand's phase is given by, 

9 = -n' u _ - u+U' - uJV' - m! j-, (69) 



and the wave packet trajectory (i.e.. dQ/dcu^ = 0), using relation (|6?|), is: 

Q + U' = in' + V) (70) 

Let us discuss these trajectories. On the "roof U' — the trajectory QT0"D is V = — n' = 
const < 0, the coordinate V is always negative, see ([T3|), thus n' is positive and, for a fixed 
wave packet, constant. This relation means that on the "roof U' = the wave packet 
trajectories are straight lines with V = const. = —n'. Similarly on the "roof V = 0, 
the trajectories, (ffO| ), are {Cj + /CjJ)U' = n', this means that they are straight lines with 
coordinate U' = const. = (uj_/uj + )n'; now the label n' is negative since the coordinate U' 
is negative, see (p!3|) . Therefore wave packets with n' positive are localized on the U' = 
"roof of the horizon and wave packets with n' negative are localized on the V = "roof 
of the horizon. These "roofs" of the horizon are depicted in Fig. 3. 



4.6 Bogoliubov coefficients 

Using the well defined inner product on the surface of the horizon ( j59f) , we can com- 
pute the Bogoliubov transformation coefficients relating the propagated "in" wave packet 
modes from the flat region (region IV) up to the horizon of region I and the "out" wave 
packet modes defined on the horizon. The two Bogoliubov coefficients, defined in the 
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USUal way (i.e.. OLi x ~ k _ nqm .^_ n , Vrnl — { v % c J t _ nqm -, u T-,n',l',m') an d ^k x ,k-,n,q,m;u>-,n',l',m' 
~( U k x ,k^,n,q,m' U C^*n' J' ,m>)) are ' 



ry~ - 1 I rk x +e rk—+8 ftjj—-\-e' 

k x ,k-,n,q,m;u,-,n>,l>,m> I _ -p_^_ / ^ / ^ _/ d(jJ _ e ~ink x -iqk-±in u- (y-Q 



@k x ,k— ,n,q,m;6j- ,n' ,1' ,m' J V'^^f' •'fc 

x i(2M) (±1)™^^ |cV]e ie «r(l + z4M|A; a; |)e ±27rM l fe *l (2Mcu ± p 4M|M 
^7r\k x \uj- 



where I' ~ y^a and 6* a are given by (|55|) and (^T]) respectively. The upper signs in (|71~1) , 
except for the label u;±, stand for the a coefficient and the lower signs for the f3 coefficient. 
The label u)± is the same for both coefficents and its double sign means k x > for the upper 
sign and k x < for the lower sign. 

In the approximation of small integration intervals e, e' and 5, and following the stationary 
phase technique one can take the non oscillating terms in the integrals as constants. For that 
reason we can perform a convenient separation in the Bogoliubov coefficients into modulus 
and phase as, 

"te-™,rn-fi,-, n >w 1 = T < 2M ) ( ±1 ) m )r(1 + iA M\k x \)\ e ±2 * M ^K a/p , (72) 

r-^k x ,k-,n,q,m;LO-,n' ,V ,m' ) \/TY\k r \ 



where we have made use of Stirling's formula for the gamma function |22 

1 1 



i.e. 



T(z) = z z ~h~ z V2^ 



+ ^r^ + 0(z 



-3\ 



12z 288z 2 



it can be written as 



r(l + iy) = |r(l + iy) I exp i (| ln(l + y 2 ) - y + i tan" 1 yj (73) 

for values y > 0, and can be written as 

T(l + iy) = \r(l + iy)\exp(y\ny-y + ir/4), (74) 
for y > 1. We have also defined 1Z a /p as, 

°rn,±m> / „ - in k x +iQ / >, \ n i iq k_+i9 a 



K a/f) = -^p^ I dk x e~ mk * +lU I dk- \C V \ e - iqk - +l " a (75) 
V e5e' Jk x 



w +e dU- ±in'ui--i4:M\k x \ln(2Mu>±) 
i ' 



where Q = 4M\k x \ In (4M|fe z |) — 4M|A; 2 .| + 7r/4, it comes from Stirling's formula Q73|), where 
we have assumed 4M|fc x | > 1 for simplicity. This assumption does not afect our results 
because the purpose is to relate a single "in" wave packet with a fix label k x to a single 
"out" wave packet with a fix label u)_. Of course we can follow the same steps for the case 
4M|fc. r | > using the more accurate version of the Stirling's formula (|73[ ). We have kept 
the terms \C\>\ and (c<j_) _1 / 2 in the integrands for convergence. 
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The terms, H a and Hp act like localizing terms in the sense that the three integrals, of 
which they consist, have a well localized peak. The integrand's phase of H a and Tip (Q a 
and respectively) are 

Q a/ p = -nk x - qk^±n'u_ + Vt + 6 a - AM\k x \\n(2Mu ± ), (76) 

and the peaks of H a and Tip are in the spacetime points which satisfy dQ a /p/dk x = 
dQ a /p/dk- = d® a /p/duj- = 0. Solving these three equations for Q a we find, respectively, 

„ = <*Lfc4Mn ( 2 M) , (77) 



dk x \ Cb± 



' a 



dk_ 



(78 



ri = -AM—, (79) 

where the upper sign in (|77|) stands for k x > and the lower sign for k x < 0. Similarly, from 
the equations for Qp, we obtain: 

n = ^±4Mln I ^ I . (SO) 



dk x \ lu± 



(81) 



n > = 4M^, (82) 

where the two first equations are functionally the same as ( |77| ) and (|78| ) but with a bar over 
the "in" labels, i.e. k x , n, q, since these three equations come from the (3 coefficient which 
relates an "out" wave packet with an "in" anti wave packet (i.e. a wave packet constructed 
by the superposition of monocromatic modes of negative frequency). On the other hand 
(O)' (HH) an d (0) come from an a coefficient and so relate an "out" wave packet with an 
"in" wave packet. In fact, notice that equations (|77|) , (ff8|) and fl79D give a relation between 
the labels of the "in" wave packets (i.e. k x , k y , k_, n and q) and the labels of the "out" wave 
packets (i.e. u_, n', I' and m!) and equations (p0|), (|81|) and (|82"D give a relation between the 
labels of the "in" anti wave packets (i.e. k x , k y , A;_, n and q) and the labels of the "out" 
wave packets (i.e. u_, n', I' and m'). When the previous equations are satisfied the terms 
lZ a and Tip have a peak, otherwise they roughly vanish. 

Let us now extract some more information from equations (|77j)-(|8"2|). For instance, equa- 
tion ([77]) (or (|8TD) can be written as, 

It I _ 1 . +xm(o)/4Af _ / -w+f/;/(4M); ^. > 0, , . 

\ k *\-r ±e ~ { -u;_y '/(4M); fc x <0, (83) 

where we have used (|5l|), (|64|) , (|65|) and the fact that n + (k x /k^)Lif(uo) = 2^(0) is the 
value of the coordinate x when the "in" wave packet coming from region II crosses the 
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boundary {v = 0, < u < vr/2} into region I (see fl3iP ). Note that U' Q or Vq are the 
coordinates on the horizon reached for the wave packets. From (|8"3| ) equations ( |77| ) and ( |HD| ) 
can be seen as redshift formulae, because they relate the energy of an "out" wave packet (i.e. 
ujJ) with the energies of an "in" wave packet and an "in" anti wave packet, respectively, 
(the energy label for the "in" wave packets in the region near the horizon is \k x \). Note 
also from ([33l), that the redshift coefficient is given by the value of the coordinate U or V 
on the horizon reached by the "in" wave or anti wave packets and by the surface gravity of 
the horizon, i.e. (4M) _1 . Equations ( ff9|) and ( j82|) are position formulae because they give 
information on the trajectories of the "in" wave and anti wave packets which are related to 
a given "out" wave packet. In fact, the label n! appearing in these formulae can be written, 
in analogy to flS3|), as 

{Co _l 
-— U Q ; k x > 0, 
uj- (84) 
-V ; k x <0. 

Recall that the label n' of the "out" wave packet is related to its position on the horizon, 



see (|70|) , when n! > the "out" wave packet is localized on the U' = "roof at coordinate 

V = V^ ax = —n' = const, and when n' < it is localized on the V = "roof at coordinate 
U' = f/^ax = (u-/Cj+)n' = const. From (|79|) we see that when nl > then k x < 0, and 
viceversa, and from fl82|) , that when n! > then fc^ > 0, and viceversa. Equation (|78| ) 
(or (|81~D ), with the use of (|5l|) , is the same as (|6lD and it does not give us any additional 
information. 

Putting all this together we can give the following interpretation. An "out" wave packet 
(a)_, n'), which is localized at the coordinate on the "roof U' = {U' max on the "roof 

V = 0) of the horizon, is related to an "in" wave packet with momentum along the x axis 
k x = o;_V^ iax /4M {k x = —Cj + U' max / AM), which reaches the horizon at coordinate Vq = V^ ax 
on the "roof U' = {U' Q = U^ ax on the "roof V = 0), and to an "in" anti wave packet 
with the same momentum along the x axis, but with opposite sign, k x = — £>-V^ ax / 'AM 
(k x = oj + U' mSkX /AM), which reaches the horizon at coordinate U' = (^-/^+)Knax on the 
"roof V = (Ift = (u>+/u-)U^ on the "roof U' = 0). 

4.7 Particle creation 

Following the formalism of quantum field theory on curved spacetime, spontaneous particle 
creation is directly related to the f3 Bogoliubov coefficient. In fact, the number of "out" 
particles in a given wave packet mode with labels (£>_, n' , V , m') (i.e. the number of quanta 
in the wave packet mode (u)_, n', I', m')) in the "in" vacuum is given by the sum over the 
"in" labels of the squared modulus of the (3 coefficient |p3|, (f72|), that is, 



(0, in|iVg nW |0, in) = _ ^ Pk x ,~k-,m,n, q -,Q-,n>,l>,rn> • ( 85 ) 

Here N~^ n , v m , is the number operator of "out" particles defined in the standard way as 

^%',V,m' = a t^n',l',rn' a l n -,n',l',m', where a ^n',l',m' and ^n'.i'.m' are tlle "° ut " WaVe Packet 

creation and annihilation operators respectively. With these operators we can write the field 
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operator 0, as a combination of "out" wave and anti wave packets, i.e. 



00) 



E 



a ^,n',l',m' U> ^,n',l',m'( X ) + °i_ ,n' ,V ,m' U ^n' ,V ,m' ( X ) 



To evaluate the sum in (85) it is worth noticing that the Tip coefficient in (72), given by 



([75]), satisfy the following two equalities: 

£ ^ 2 = ww (86) 

£ l*f>| 2 = -^rjvv ( 87 > 



iD_ .n'l'm' 



k x ,k—,m,n,q 



4(2Mf 



These equalities follow from the general property (§) of wave packets and from the form of 
Cy given by flBCp. Note also, in order to make sense of the square of Cy according to ( JODj ) 
we consider the product of a Dirac's delta with a Kronecker's delta. 

These two equalities and the fact that the Tip term has a peak when the relations (|80|)- 
([32]), between the "in" labels (i.e. k x , k y , n, g) and the "out" labels (i.e. n', i', m') 
are satisfied, allow us to approximately write 

^ ~ 2(2M) 3 / 2 ^."'"-C^) -1 ^n,n(u>-,n',l',m') $k-,k-(u>-,n',l',m') $q,q(u>-,n',l',m') ^m,-m'i (88) 

where n = n(uj-,n',l',m') is given by (]8"0|), g = q(u}_,n',l',mf) by (]8~1[), and „/,j', m ') 
is given by ^' j / a which appears from the squared modulus of (|60|). Then the number of "out" 
particles created in the "out" packet mode (uj-,n',l',m'), equation ( [Sip , is simply, 

(0, inl^y^^lO, in) = N^ n , = e8wMQ _ (]n y 4M) _ ^ (89) 

which can be interpreted as a thermal spectrum for each fixed value of the label n', with a 
temperature, 

T = r. (90) 

8ttM(H/4M) v ; 

This spectrum depends on the dimensionless label n'/4M, i.e. on the trajectory of the wave 
packet, but it is thermal for all wave packet modes with the same trajectory. This is quite 



different from the black hole case [0] where the temperature is independent of the packet 
trajectory ]17j and depends only on the surface gravity, k = (4M)™ 1 . In ref. [[J we discuss 
how the black hole case can be seen in some sense as the time reversal of the colliding wave 
case. The physical interpretation of this n'-dependent temperature follows from the fact that 
we are computing the particles produced on the wave packet mode (u)_, n', I', m'), which is 
localized on the horizon by n', given by ([TO]), so that these particles may be "localized" in 
the same position on the horizon. Note that when n' ~ 0, i.e. near the bifurcation point 
U' — V — of the horizon, the temperature is higher. 

This spectrum is in agreement with the spectrum of particles created on the monocro- 
matic modes (|20D, i.e. particles with a well defined momentum but not localized in space, 



because the one particle Fock space can be decomposed on a basis given in terms of the 
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monocromatic labels, \u), or, alternatively, on a basis given in terms of the wave packet 
labels, \oj,n'). The space is the same but, of course, the particle interpretation is different. 
For the monocromatic modes, discused in ref. ||, the spectrum of particles created (|20|) is 
inversely proportional to the energy of the "out" modes (i.e. u;_), with a proportionality 
factor (87rM) _1 . The relation between the number operators and N^ n , is given in 
terms of the Bogoliubov coefficients relating the monocromatic modes and the wave packet 
modes. Such transformation has no (3 Bogoliubov coefficient because the positive frequency 
wave packets are constructed with positive frequency monocromatic modes only. The a 
Bogoliubov coefficient is 



a 



ui- \lo— ,n' 



where, here, 9{x) is the usual Heaviside step function, thus the relation between the monocro- 
matic and wave packet annihilation operators is given by 



11 s — i(n'— n')u>- „t out n out 



Then the monocromatic number operator can be written as 

and we can approximate 

(0, in|iVr|0, in) ~ I X)(0, H«R<?,J0, ™> = ^ £A&, 

n' n' 

as the phase term exp[— i(n' — n')uiJ\ is fast oscillating (except when ml = n'), because 
n' — n' = (I — l)2ir/e' with I and I integers and e' a small positive parameter. Then 

9M 1 roo Ah 
(0, in|iVr|0, in) * ^ -4— \ (91) 



7i 8ttMuj- Jo e 

where we have assumed 8ttMuj_ « 1 in order to approximate k = (&7iMujJ)n' j 'AM as a 
continuous dimensionless variable. The inverse proportionality in uo_ and the logarithmic 
divergence, in fl2~0|), have been recovered. The extra factor (2M) 2 in ( p0|) is due to the fact 
that two monocromatic labels in ref. || were, in fact, discrete and this factor guarantees 
the correct normalizations. 
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Figure Captions 



Fig.l 

This figure shows a projection of the colliding plane wave spacetime in the (u, v) plane. 
One can see four regions: region IV is the flat region before the waves collide, region II and 
III are the plane wave regions and region I is the interaction region. The interaction starts 
at u = v = 0, the lines u = 0, v < and v = 0, u < are the boundaries of region IV with 
regions II and III respectively, and v = 0, < u < n/2 and u = 0, < v < n/2 are the 
boundaries of region I with regions II and III respectively. Atu + t> = 7r/2a Killing- Cauchy 
horizon is formed which one can see as a coordinate singularity of the metric in region I. 

Fig.2 

In this figure we represent the coordinates (£, x) in terms of the Kruskal-Szekeres like 
coordinates (£/', V) in the interaction region. The lines £ = constant are hyperbolae and 
the x = constant are straight lines crossing the origin U' = V — 0. The Cauchy horizon is 
{V = 0, V < 0} U {V — 0, U' < 0} which corresponds to the limit of the hyperbolae as 
£ — > 0. The "roof U' = corresponds to x — > — oo and the "roof V = to x — > oo. 

Fig.3 

This is a 3-dimensional plot of the interaction region (region I) in which all its boundaries 
are shown, using non singular Kruskal-Szekeres like coordinates. The surface S3 and the lines 
.Mi, M.2, -M'i and M.' 2) are the boundaries of region I and the two plane wave regions II 
and III, the points V and V' correspond to U' — 0, V — 0, r\ — —tt/2 and to U' — 0, V — 
0, rj = 7r/2, respectively, and they are folding singularities P, 10, ^4|. The Cauchy horizon 



is the "roof {U 1 = 0, V < 0} U {V — 0, U' < 0} and region I is enclosed between the 
surface S3 and the "roof. 
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